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^ ABSTRACT 

We investigate the effects of a top-heavy stellar initial mass function on the reionisation his- 
tory of the intergalactic medium (IGM). We use cosmological simulations that include self- 
consistently the feedback from ionising radiation, H2 dissociating radiation and supernova 
(SN) explosions. We run a set of simulations to check the numerical convergence and the 
, effect of mechanical energy input from SNe. In agreement with other studies we find that it 

^ ■ is difficult to reionise the IGM at Zrci > 10 with stellar sources even after making extreme 

assumptions. If star formation in 10^ galaxies is not suppressed by SN explosions, the op- 
i tical depth to Thomson scattering is ^ 0.13. If we allow for the normal energy input from 

CO ' SNe or if pair-instability SNe are dominant, we find Tc ^ 0.09. Assuming normal yields for 

I the first stars (population III), the mean metallicity of the IGM is already Z/Zq — 2 x 10^^ 

(10^'^ < Z/Zq < 1 in overdense regions) when the IGM mean ionisation fraction is less than 
10 %. For these reasons population III stars cannot contribute significantly to reionisation un- 
less the mechanical energy input from SNe is greatly reduced and either the metal yield or 
• the mixing efficiency is reduced by a factor of 10''. Both problems have a solution if popu- 

Qh| lation III stars collapse to black holes. This can happen if, having masses A/* < 130 Mq, 

they are characterised by heavy element fall -back or if, having masses > 260 Mq, they 
^ ■ collapse directly onto black holes without exploding as SNe. If metal-poor stars are initially 

^ I important and collapse to black holes is the typical outcome, then the secondary emission of 

. ionising radiation from accretion on SN induced seed black holes, might be more important 

than the primary emission. 

We also develop a semianalytic code to study how Tc is sensitive to cosmological parame- 
■ ters finding essentially the same results. Neglecting feedback effects, we find simple relation- 

5-H ' ships for Tc as a function of the power spectrum spectral index and the emission efficiency of 

ionising radiation for cold dark matter and warm dark matter cosmologies. Surprisingly, we 
estimate that a warm dark matter scenario (with particle mass of 1 .25 keV) reduces Tc by only 
approximately 10%. 

Key words: cosmology: theory, dark matter - galaxies: dwarf, formation, haloes - methods: 
numerical 
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1 INTRODUCTION 

It is widely believed that the earliest generations of very low metal- 
licity stars would have been far more effective than normal stars 
in reionising the universe, so much so that they could account for 
the large optical depth to Thompso n scattering, Te = 0.17 ± 0.04, 
measured by the WMAP satellite ( Bennett et alllOOBllKogut et alj 
[2003). The main reason for this expectation is that the stellar initial 
mass function (IMF) is thought to be tipped towards high masses 
for Z/Zq <g 10"*. But there are two byproducts which are likely 
to follow from this scenario. The same stars, which are efficient UV 
producers, will die as supemovae (SNe), the explosions heating the 
surroundings and inhibiting further star formation ("negative feed- 
back"), and the metals ejected in these explosions will contaminate 



the high density regions, rapidly bringing the metallicity up to a 
level where such stars cannot form. We argue that, while it is pos- 
sible to evade these strictures, it is very difficult to do so unless this 
early generation ends its life primarily via implosion to black holes 
rather than explosion as supemovae, and that in the former case 
the secondary effects of ionising radiation from accretion onto the 
formed seed black holes may dominate over the primary UV from 
the first generation stars. 

In this discussion we are focusing on the period 15 < z < 25 
when the haloes within which star fo rmation would occur are 
small and gas cooling is initially via H2 iCouchman & ReeJl986t 
Ostriker & Gnedin 1996). The later phase, responsible for the more 
complete reionisation of the universe as attested to by the Sloan 
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Digital Sky Survey (SDSS) quasars at 6 < z < 7 iBecker et all 
l200lt lFan et al. 2003), is fairly well understood and can easily be 
accounted for by a relatively normal population II. For a Salpeter 
mass function of population II stars in the mass range IM© ^ 
M* ^ IOOM0 , and an escape fraction of UV p hotons from th e 
forming galaxies to the IGM of (/esc) = 20 - 50%.'GnedirJ j2000l) : 
[Ciardi et al. (2003); Ricotti (2003a); Chiu et al. (2003) among oth- 
ers, find little difficulty in matching the QSO-determined ionisa- 
tion history using standard stellar interiors theory and cosmologi- 
cal models. Th e two epochs of ionisation might have been distinct, 
as proposed bv lCenl i2003 hf). or there may have been an extended 
period of partial ionisation with a gradual transition to the epoch 
observed in quasars. If the IGM had a sudden transition from neu- 
tral to completely ionised, the measured optical depth To = 0.17 
corresponds to a redshift of reionisation Zrci = 17 for the best fit 
ACD M cosmological model measured by WMAP iSoergel et alj 
l2003h . 

The large Tc measured by WMAP implies that an additional 
source of ionising radiation is important at high-redshift . Two 
categories of sources have been considered in previous studies: 
(i) stellar or AGN sourc es, forming at z 40 in the first galax- 
ies (e g.. .Cea ■2003a; .Wvithe & Loebl 120031; ISomerville & Liviol 
l2003t ICiardi et alj|2003l ; ISokasian et alj|2003h or (ii) exotic (un- 
known) sources such as decaying neutrinos (e.g.. ISciamalll983 ; 
[Hansen & Haiman 2004), super-heavy dark matter decay at recom- 
bination (Bean et al. 2003) or primordial B Hs, either via evapora- 
tion ( .Guedens et aU2002) or via accretion jGnedin et alfiggS) . In 
this paper (paper I) we study the early reionisation of the IGM by 
the first stars (population III). In two companion papers (paper Ila 
and p aper IIb,fRicotti & Ostriker 2003; Ricotti. Ostriker. & Gnedin 
l2004h we study the partial ionisation by a putative X-ray back- 
ground produced by mass accretion on seed BHs formed in the first 
galaxies. 

Several g roups have studied, using analytical calculations 
(e.g., lUehara et al. 199 6 ; Larson 199& ; .Nakamura & U memura 
I999h or numer ical simulations {e~, lAbel et al] 120021 ; 
Bromm et alJ Il999h . star formation in a metal-free gas. Most 
groups agree (but see Omukai & Palla 2003, for a more conserva- 
tive point of view) that the first population of stars was probably 
massive, with typical mass of about 100 M©. The stellar IMF in 
the early universe is unknown but it is e xpected to be top-heavy or 
bimodal iNakamura & Umemural200ll) . It seems natural therefore 
that at high-redshift the emissivity of ionising radiation was larger 
than euv ~ 1.3 x 10~* found for a Salpeter IMF. If all stars 
were massive (M, ^ 10 M©), the emission efficiency of ionising 
radiation would be a factor of 10-20 larger than for the Salpeter 
IMF. During their lifetime massive stars convert a little less than 
half of their hydrogen into helium. The efficiency of thermonuclear 
reactions for the p-p chain is about 7 Mev per hydrogen atom. 
If we assume that all the energy liberated is emitted in ionising 
photons we obtain a theoretical upper li mit for the UV emissio n 
efficiency e'Q^ < 2 - 3 x 10"^ (c/, iMiralda-Escude "20031). 
The maximum emission efficiency eyv" can be obtained only if 
low-mass stars do not form and if all the ionising radiation can 
escape from the molecular clouds and the diffuse ISM of the host 
galaxy. Here euv is the ratio of energy density of the ionising 
radiation field to the gas rest-mass energy density converted into 
stars iptC?). Thus, euv ~ {hpV /mH(?)N^^, where A^uv is the 
emitted number of ionising photons (with mean energy hpv) per 
baryon converted into stars. At first glance, the enhanced efficiency 
of UV production from population III stars seems to be the natural 
explanation for the WMAP result. But as noted earlier there are 



several issues that need to be addressed that make such a model 
difficult to achieve; 



• Complete reionisation at increasingly high redshift requires a 
larger number of ionising photons per baryon. The recombination 
rate of hydrogen is large and the fraction of radiation escaping each 
galaxy is likely to be less than unity. 

• The efficiency of star formation in small-mass haloes (A/dm < 
10** M0) is regulated by radiative feedback effects (i.e., star for- 
mation is reduced as a result of gas heating by UV, winds and 
SNe produced by star formation) and, according to cosmological 
simulations (Ricotti et al. 2002a b), these galaxies cannot reionise 
the IGM. The simulations have shown that even assuming a top- 
heavy IMF, the star formation is strongly suppressed and the mean 
stellar mass in small-mass galaxies decreases. The simulation re- 
sults agree with observations of the smaller-mass galaxies (dwarf 
spheroidal galaxies) in the Local G roup that show the same large 
mass-to-light ratios <Ricottil2003bh . 

• Later forming galaxies with masses Mdm ^ 10^ Mq have a 
larger star formation efficiency and their stars could reionise the 
IGM. But the enhanced energy injection by SN explosions due to 
a top-heavy IMF would also reduce star formation in larger galax- 
ies. If the energy of the SNe is normal (£51 = -E/[10'^^erg] = 1), 
the SN energy input for a top-heavy IMF would be larger than for 
a Salpeter IMF. Some massive stars are known to explode as hy- 
pernovae with energies E^i ^ 1. Supermassive stars with masses 
I4OM0 < < 260 Mq are thought to end their lives with 
an explosion (pair-instability SN) that does not leave any remnant. 
The explosion energy of pair-instability SNe (PIS) is proportional 
to their mass and is typically ^ 100. 

• SN explosions eject metals into the ISM and IGM. The tran- 
sition from a top-heavy IMF to a Salpeter IMF is thought to de- 
pend on the ^as metallicity. If the metall icity is Z ^^7 ^ ~ 10~^ 
Zq tomuk ai 2000*; 'B romm et al .''2001; Schneider et al. 2002) the 
gas can cool more efficiently and to lower temperatures, allowing 
fragmentation and the formation of small-mass stars. The impor- 
tance of population III stars for reionisation depends on how long 
metal-free gas is present in protogalaxies. This issue is difficult to 
address because it requires a good understanding of the mixing and 
transport processes of the heavy elements, but metals produced by 
SNe will clearly lead to an enviro nment in which a top-heavy mass 
function no longer exists (e.g., see IWada & Venkatesanl2003 h. The 
only ways to avoid contamination that we have thought of are two: 
(i) have such efficient ejection of metals from star forming regions 
that lower density regions are contaminated before higher density 
ones (d\n Z / Zq / dln(p) < 0) or (ii) have such inefficient local 
mixing that effective homogenisation is delayed for 0.5 — 1 Gyrs, 
~ 5 — 10 local dynamical times. We will return to these possibilities 
later. 



This paper is organised as follows. In §|2|we show, using order 
of magnitude estimates, that a large To favours a scenario in which 
a large fraction of the first massive stars collapse into primordial 
massive BHs. In §|3|we show the results of a set of cosmological 
simulations with radiative transfer We focus on the effects of the 
limited mass resolution of the simulations and feedback from SN 
explosions. In § 13 we use a semianalytic calculation to estimate 
the effect of cosmological parameters on To. In §|5|we present our 
conclusions. 
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2 BH-FORMING OR PAIR-INSTABILITY 
SUPERNOVAE? 

In this section we use simple arguments to estimate the number of 
ionising photons needed to reionise the IGM. Since, conventionally, 
the massive stars that produce ionising radiation end their lives as 
SN explosions and BHs, we estimate the number density of metals 
and relic BHs expected given the redshift of reionisation (or To). We 
show how the observed metallicity of the IGM and ISM in galaxies 
and the demographics of supermassive black holes (SMBH) at 2 = 
can put some loose constraints on the IMF of the first stars. We 
find that it is very difficult to reionise the universe at high redshift 
without strongly contaminating it with metals. These calculations 
admittedly rely on simplistic assumptions but are useful to show 
the problems related with each scenario. 

Assuming ionisation equilibrium we have, 

dt t\l {tree) 

where riph is the number density of ionising photons, Uc is the 
mean electron number density, tu is the Hubble time and {t^cc) ~ 
(aBTieCHii)'^ is the mean recombination time inside the Hll re- 
gions. The recombination and the Hubble times are, respectively. 



(1) 



{tree ) 



1.7Gyr 



IGyr 



0.02 



7 



1 + z 
7 

h_ 
07 



^-1 



Te — car 



dto, - 0.069^^ J dz{l + z)'-^^{xe)M, 



ph 



nn {t icc) car 



10 



0.1 



ph 



H II regions it is not constant with redshift because it depends on the 
spatial distribution of the sources that are responsible for reionisa- 
tion. The effective clumping factor, Chu, is larger for reionisation 
at higher redshift because the sources, being more numerous and 
with smaller mean luminosity, must ionise first a large fraction of 
the dense filaments in which they are located. In the following cal- 
culations we use equation j6j because we want to use the most con- 
servative assumption on the number of ionising photons A^ph. Us- 
ing equation Q would only change the exponent of Tc from unity 
to four in the following relationships of this section. 

If the massive stars that produce ionising radiation end their 
lives exploding as SNe and eject most of their mass enriched with 
heavy elements, the metal yield, Y, of a star population is approx- 
imatively proportional to the number of ionising phot ons em itted, 
independently of the assumed IMF (e.g., Madau & Sh ulll99d). Us- 
ing the population synthesis model STARBURST ( L eitherer et alJ 
11999), we found that a population of stars with metallicity 
Z ^ 0.05 Zq (population II) has a yield 



(2) 



(3) 



where Chu = (whii)/(whii)^ is the mean effective clumping fac- 
tor of the ionised gas. The number of ionising photons emitted per 
hydrogen atom (also counting the number of recombinations per 
Hubble time) sufficient to fully reionise the IGM at redshift z^ai is 
given approximately by 



ph = ~ ma.x{tH /{Uoc), 1) ~ max([(l+Zrci)/14]^'"'''CHii, 1), 

(4) 



where nn is the mean hydrogen density. Assuming sudden reioni- 
sation at redshift Zrd > (UC^^/^ — 1) and neglecting the redshift 
dependence of the clumping factor, we have A'^ph oc (1 + ZrciY'^- 
The optical depth to Thomson scattering is given by 



(5) 



where {xs)m is the mass weighted electron fraction. Assuming a 
sudden reionisation at redshift Zrci we have Tc oc (1 + Zrei)^'^, that 
has the same redshift dependence of A'^ph. Therefore, in the case 
we neglect the redshift dependence of the clumping factor in equa- 
tion J4j, we find t hat A^ph is proportional to Tc. A similar relation- 
ship was found bvloh et alj i2003) that, owing to the cancellation 
of the electron density in equation jsj, also find the following ex- 
pression: 



(6) 



that it is independent of redshift. The results of numerical simula- 
tions, presented in S I3.2I will show that this expression is not very 
accurate. A better agreement can be obtained using the equation 



(7) 



that has been empirically derived from the simulations. The reason 
for this result is that the mean effective clumping factor inside the 



Y = 155 Zg 



2 X 10- 
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where 0.3 < g < 2 depends weakly on the metallicity and euv is 
the efficiency of ionising radiation emission from stars. As noted 
in the introduction, the maximum efficiency for thermonuclear re- 
actions is euv" ~ -Bp-p/mnc^ ^2 - 3 x 10"^, where -Bp-p is 
the mean energy per nucleon emitted in form of ionising radiation. 
In this calculation we have assumed that population III stars during 
their lifetime are able to convert half of their hydrogen into helium 
and that they are hot enough to emit all their energy into ionis- 
ing photons (cf., Miralda-Escude 2003). The yield of population III 
is very uncertain. The main uncertainty comes from the unknown 
fraction of their stellar mass and metals that implodes into a black 
hole. In the following paragraph we will show that, if population III 
stars have "normal" yields (neglecting that a substantial fraction of 
their mass and metals is locked into the black hole remnant) they 
cannot produce the large Tc measured by WMAP without a produc- 
tion of metals that is not consistent with the observed metallicities 
in the IGM and ISM. 

The emissivity of ionising photons by massive stars is 



A:(1.36 X 10"**)Aru, 



whereX — hu /13.6eV ~ 2 — 3 is the mean energy of H I ionising 
photons and Nyy is the number of ionising photons per baryon 
converted into stars. If /, is the fraction of baryons converted into 
stars we have A^ph = N^^f*- 

Let's assume that massive population III stars have the same 
yields as massive population II stars. Since the total metal produc- 
tion is Z = Y ft, using equation Js}, we have 



6.8 X 10"'' i ^ ^ ^° I XA^ph and 
euv 

{Ig X 10-" Zo)XA^ph. 



(9) 
(10) 



The fraction of baryons that needs to be converted into stars to get 
a given Tc can be estimated using equation J6} and using X = 3: 



/, ~ 2 X 10" 



2 X 10"^ 



0.1 



This relationship shows that, if the IMF is top-heavy, a smaller frac- 
tion of baryons needs to collapse into stars. We can also estimate 
the mean star fraction of galaxies /f^' — /. //coii, where /coii is 
the fraction of baryons in virialised DM haloes. At 2 = 16, most 
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of the collapsed baryons are in small-mass galaxies. The mass frac- 
tion of DM haloes with mass Mdm < 10* M© is fidm ~ 1-74 % 
and, fldm ~ 0.05 % in haloes with mass Afdm > 10* M©. If we 
assume that the baryons follow the DM (/coii = fidm), the fraction 
of gas converted into stars must be 




(^) if A/dn. > 10*Mo 

{ ^) if A/dm < 10* Mq. 



(11) 



If galaxies in small-mass haloes do not contribute to reionisation, 
the efficiency of star formation in larger galaxies needs to be ap- 
proximately 30 — 50 %, which is perhaps larger than the efficiency 
estimated for globular clusters. If small-mass galaxies contributed 
significantly to reionisation, then their efficiency of gas conversion 
into stars must be larger than 1 — 2 %. Assuming a "closed box" 
chemical evolution model for the galaxy, the ISM metallicity is 



^isM = F/r' = 



6»Zo(^) if Mdm > 10* Mo 
0.155 Z0 (^) if A/dm < 10* Mg 



(12) 

These metallicities are far too large when compared to the metal- 
licities of the oldest star clusters known: old globular clusters that 
typically have [m/H]=-I.5 {e.g., Vandenbera et al. 1996) and dwarf 
spheroidals that have -2.5 < [m/H] < -1 (e.g.. iMatedll998i) . To 
match observations all but 10~^ to 10~^ of the metals would need 
to have been ejected. With such a large in situ metal production it is 
also difficult to have patches of metal-free gas so that population III 
stars can continue to form. If instead we assume that most metals 
are ejected from the galaxy haloes into the IGM, we find (assuming 
X = 3) that the IGM mean metallicity would be 



ZiGM = 3g X 10 Zg 



0.1 



(13) 



But this value is about ten times larger than the estimated mean 
metallicity Z ~ 2 x 1 0~^ Z© of mildly o verdense regions in the 
IGM at redshift 2: ~ 3 jSchave et all2003h . It seems therefore that 
we need to find a mechanism to hide most heavy elements. If metal 
enriched gas remains highly inhomogeneous in the ISM, popula- 
tion III stars could continue to form for a while, together with nor- 
mal stars. This scenario is not implausible but, because of the for- 
mation of some small-mass stars, the efficiency of UV emission 
will be reduced to euv < 10"'' and reionisation would not occur 
by this process. 

Let's now relax the assumption that population III stars have 
the same yields as population II stars. The initial mass function and 
metal yields from population III stars are extremely uncertain^ The 
simplest way to parametrise these uncertainties is to assume that a 
fraction, /bh, of population III star mass and metals collapses onto 
BHs. The baryon fraction in seed BHs, cijbh = /bh/*, and the 



^ . Metal-free stars in the mass range IOMq < Af, < 130 have 
reduced yields because, being more compact than population II stars, they 
are characterised by a large amount of heavy element fall-back onto the 
black hole remnant. Stars more massive than 260 do not explode as 
SNe but rather collapse directly into a BH without exploding. On the other 
hand, stars with masses 13OM0 < M» < 260 M0 provide the worst 
case for metal pollution and mechanical energy input. They explode as pair- 
instability SNe (i?5i ~ 10 — 100), ejecting all of the metals without leav- 
ing any remnant and leaving a specific signature in the relative abund ance 
distrib utions that has been sought for but not found (.Venkatesan & Truranl 
l2003ft. 



metal production are given by 

'2 X 10"^ 



<J-'BH 



/bh2 X 10- 
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Z ~ (1 - /bh)33 X 10-^ Zq 



(14) 



(15) 



An upper limit for the total mass of seed BHs can be estimated 
assuming that most of the mass in supermassive black holes 
(SMBHs) observed in the bulges of galaxies at 2; = is built 
up from the merger of seed BHs. We estimate a baryon frac- 
tion in SM BHs, tJBH.max, from the BH-bulge mass relation {e.g., 
iKormendv" & Richstone 1995) and from the estimated baryon frac- 
tion in spheroids, ojbuigc, at 2 = 0: 

/ A'fsMBH \ ,1 I ^^ r,\ 1^-4 

' ' t^buigo = (1 ± 0.7) X 10 



WBH,max 



\ A'/bulgc / 



We adopted the value s A^sMBn/Afbuigo = (1. 5 ± 0.3) x 10~^ 
jGebhardt et al.""2000) and (.Jbuige = (6.5 ± 3)% (Per sic & Saluccil 
[1992 ; Fukuaita et al. 1998). This is an upper limit because a sub- 
stantial fraction of the mass of SMBHs must be due to gas ac- 
cretion to account for the luminosity of quasars. It is also possi- 
ble, however, that a fraction of seed BHs is expelled from their 
host galaxies or does not end up in the observed population of 
SMBHs iMad au & Rees 2001..) . The recent discovery of several 
ultra-luminous X-ray sources (ULX) has been interpreted as evi- 
dence for intermediate-mass BHs (of about 100 Mq). These objects 
could be relics of primordial BHs. Theoretical studies have shown 
that the expected luminosity function of accreting r elic BHs in 
galaxies is consistent with observations of ULX source ilslam et all 
2003a b). Since we do not know how much they contribute to the 
BH mass budget at z = in the calculation presented here, we 
neglect their contribution. In any case, in paper II and paper III we 
will show that the upper limit on usyi adopted here is valid inde- 
pendently of the fate of the first BHs because it is constrained by 
limits on the observed 7-ray background. If we use the baryon frac- 
tion in SMBHs at 2 = 0, t^BH.max = 10"'', as an upper limit for 
ui-BYi in equation <14> and Zigm = 2 x 10"^ Zq as an upper limit 
for the IGM metallicity given by Ziqm ~ {fci)Z, where (/cj) is 
the fraction of metals ejected into the IGM and Z is from equa- 
tion <15> . we find: 



1 



0.07 



T^) < /bh < 0.5 



euv 



0.1 



(7(/ej) VO.17 - ^ V2 X 10-3^ 

We estimate that, if population III stars reionised the IGM to the 
value measured by WMAP (tc ~ 0.17), the mass fraction of pop- 
ulation III stars that might collapse into BHs is /bh ~ 30 % and 
the metal yield of population III stars is less than half the "normal" 
value. This result is rather extreme but plausible if population III 
stars have masses Af, > 260 Mq or if their explosion energy is 
£51 < 1. Using the same assumptions as before on the metal and 
BH budget, we find an upper limit on the Thomson optical depth 

0.007 



Te < 0.05 



^uv 
2 X 10-3 



In conclusion, if the optical Tc ~ 0.17 of the IGM is due 
to ionising radiation from massive stars, we need to assume that 
most of the heavy elements synthesised in the stars are not ejected 
into either the ISM or IGM. Similar conclusions have been found 
in order to prevent IGM over-enrichment with heavy elements by 
population III stars, if they have to account for the observed near- 
infrared background excess after subtraction of 'normal' galax- 
ies contribution jB ond et al. 1986; Sh chekinov 1986: Santos et all 
l2002LISalvaterra & Ferraral2003h . 
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Table 1. List of simulations with radiative transfer. 



# 


RUN 


A^box 


^box 


Mass Res. 


Res. (com.) 


Smaller galaxy 


e* 


<:uv(/esc) 


FlMF 






h-^ Mpc 




pc 


h-^ Mq 










1 


64L1VM 


64 


1.0 


3.15 X 10^ 


976 


3 X 10' 


0.1 


10-3 


17 


2 


64L2VM 


64 


2.0 


2.52 X 10"^ 


1952 


2.5 X 10* 


0.1 


10-3 


17 


3 


128L1VM 


128 


1.0 


3.94 X lO'' 


488 


4 X lO'^ 


0.1 


10-3 


17 


4 


128L2VM 


128 


2.0 


3.15 X 10^ 


976 


3 X 10^ 


0.1 


10-3 


17 


5 


128L4VM 


128 


4.0 


2.52 X 10"^ 


1952 


2.5 X 10* 


0.1 


10-3 


17 


6 


128L2VMb 


128 


1.0 


3.15 X 10^ 


976 


3 X lO'' 


0.1 


10-3 


1 



Pai'ameter description. Numerical parameters: N^^^ is the number of grid cells and Lbox is the box size in comoving h^^ Mpc. The mass of the DM 
particles is listed in column four, the comoving spatial resolution in column five and the minimum DM mass of resolved galaxies (with at least 100 DM 
particles) is given in column six. Physical parameters: e» is the star formation efficiency, euv is the ratio of energy density of the ionising radiation field to 
the gas rest-mass energy density converted into stars (depends on the IMF), (/esc) is the escape fraction of ionising photons from the resolution element and 
i^lMF is a weighting parameter for the energy input and metal yields from SN explosions, that differs from unity if the IMF is not a standard Salpeter IMF 
(see S ll2l . 




30 25 20 15 10 30 25 20 15 10 30 25 20 15 10 
Z-Fl Z+1 Z+1 



Figure 1. (top panels) Star formation rate as a function of redshift for the simulations in tableQ The thick lines refer to simulations with Af^ox = 128 and the 
thin lines refer to simulations with N^^y^ = 64. All the simulations except 128L2VMb (thick dashed line) include strong feedback from SN explosions and 
differ only in box size. The Lbox = 1 h— ^ Mpc simulation resolves small-mass galaxies in which star formation is strongly suppressed by SN explosions. We 
achieve convergence in both Lbox = 2 and Lbox = 4 h"^ Mpc simulations, where the galaxies that dominate the SFR at 2 > 10 are numerically resolved. 
In simulation 128L2VMb we reduced the energy input from SN explosions by a factor of ten. The resulting SFR is about 20 times larger Integrated number 
of ionising photons produced per free electron, N-p^/ {xc), as a function of redshift (bottom panels). The redshifts after reionisation are marked with points. 
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A very inefficient mixing of the ejecta with primordial gas is 
also a possible scenario. This would lead to the situation in which 
the high density regions, preferred for subsequent star formation, 
are contaminated by metals by more than the average amount but 
the metals are not incorporated into subsequent stars formed lo- 
cally until after some delay. If the mechanism of metal ejection 
leads to metal exclusion from the star-forming high density regions, 
multiple generations of population III stars could form in the same 
galaxy and in neighbouring galaxies. In this scenario the gas metal- 
licity and overdensity are anticorrelated. Numerical simulations fail 
to reproduce this anticorrelation but this could be a consequence 
of the insufficient numerical resolution and the badly understood 
physics of metal transport. As a consequence most metals would 
reside in the Lyman-a forest at high redshift. Observation of the 
metallicity in the Lyman-a forest and the metallicity in quasar host 
galaxies at « ~ 3 do not show anticorrelation of density and metal- 
licity, but this scenario cannot at present be ruled out at high red- 
shift. 

Perhaps a more natural choice is to assume that most of the ex- 
cess metals are locked inside intermediate-mass BHs. The recently 
discovered lower metallicity star HE0107-5240 (with [Fe/H] = 
—5.3) is iron deficient and carbon rich (Christlieb et al. 2002). The 
abundance pattern of HE0107-5240 can be explained by assuming 
that the gas has been enriched by a subluminous supernova explo- 
sion (i?5i ^ 0.3) of a ze ro-metallicity star with a mass of 25 
iUmeda & Nomoto 2003). The iron metal yield is reduced because 
of fall-back into the central BH. The abundance pattern of HE0I07- 
5240 is not compatible with yields from pair-instability SNe. To 
summarise, even with extremely optimistic assumptions concern- 
ing the UV efficiency of low metallicity stars we find that the only 
scenario capable of ionising the universe at high redshift is one in 
which a large fraction of the mass of population III stars collapses 
onto massive black holes. 



a Salpeter IMF, (/esc) must be in the range 2 0-50% in order to have 
reionisation at z^d — 7 (e. e. . iBenson et air2 00 1 : Gnedin 2000). 

In this paper, as we are interested in determining an upper 
limit on Tc, we adopt the maximum emissivity of ionising radiation 
eijv' from stellar sources produced by thermonuclear reactions and 
(/esc) = 0.5. Assuming a constant emissivity is obviously not self- 
consistent since the metals ejected by population III stars eventu- 
ally trigger the transition to the era when population II stars domi- 
nate the global star formation rate. In our simulations population III 
stars are dominant for long enough to fully reionise the IGM, but 
it should be kept in mind that, as noted above, due to metal enrich- 
ment, a complete reionisation by population III stars might never 
be achieved. In the next subsection we describe the main features 
of the code. 

3.1 The Code 

The simulations were performed with the "Softened Lagr angian 
Hydrodynamics" (SLH-P'^M) code described in detail in GnediiJ 
( 1995). The cosmological simulation evolves coUisionless DM par- 
ticles, gas, "star-particles", and the radiation field in four fre- 
quency bins: optically thin radiation. Hi , Hex and Hell ionising 
radiation fields. The radiative transfer is treated self-consistently 
(i.e., coupled with the g as dynam ics and star formation) using the 
OTVET approximation iGnedin & Abel 2001.) . The star particles 
are formed when the dense cooling gas in each resolution element 
sinks below the numerical resolution of the code. The code adopts a 
deformable mesh to achieve higher resolution in the dense filaments 
of the large scale structure. We solve the line radiative transfer in 
the H2 Lyman- Werner bands for the background radiation. We in- 
clude secondary ionisation of H and He, heating by Lya scattering, 
detailed H2 chemistry and cooling, and self-consistent stellar en- 
ergy distribution of the sources cRicotti et al..2002a) . 



3 COSMOLOGICAL SIMULATIONS 

The results shown in the next section are based on numerical cos- 
mological simulations that include radiative transfer and feedback 
from SN explosions. The code has been implemented to study 
stellar reionisation and the formation of the first galaxies that is 
self-regulated by radiative feedback effects in the early universe. 
The simulations were run on COSMOS, an SGI origins 38000 in 
DAMPT, University of Cambridge. 

We adopt a concordance ACDM cosmological model with pa- 
rameters: Qo = 0.3, r^A = 0.7, h = 0.7 and f2b = 0.04. The 
initial spectrum of perturbations has as ~ 0.91 and n — 1. 

The effective emissivity is — (/csc)euv. The value of 
euv depends on the initial mass function (IMF) and the metallic- 
ity of the stellar population. Assuming a Salpeter IMF with star 
masses between IM0 ^ M* ^ IOOM0, we have for popu- 
lation III (population II) stars euv = 3 x 10~'^(1.3 x 10 '^) , 
X = h^/{l3.6 eV) = 2.47(1.76) (e.i?.. lRicotti et al.ll2002ii) . 
The quantity (/esc) here is defined as the mean fraction of ionis- 
ing radiation that escapes from the resolution element^ . Assuming 

^ Note that the definition of (/esc ) can create confusion. In numerical sim- 
ulation, such as the ones presented in this work, (/esc) is usually defined as 
subgrid physics, but in other works is defined as the mean fraction of ion- 
ising radiation escaping the galaxy haloes (e.g., the virial radii). The defini- 
tion of (/esc) and the definition of the clumping of the IGM are related and 
depend on how the outer edges of galaxies are defined. It is therefore dif- 



3.1.1 Feedback from SN explosions 

We include the ef fects of SN explosions using the method discussed 
in lGnediJ ^99^. The current resolution of cosmological simula- 
tions does not allow numerical resolution of radiative shock fronts 
and the modelling of SN explosions from "first principles". A semi- 
analytic recipe has to be included at the sub-grid level. This treat- 
ment introduces new free parameters in the simulation and the re- 
sults might depend on the model implemented. The approximation 
used to solve radiative transfer is instead much more reliable and 
has been shown, for simple test cases, to reproduce the exact so- 
lution. This is part of the reaso n why in our previou s study on the 
formation of the first galaxies <Ricotti et alj|2002d) the effect of 
SN explosions was not included. In the present paper we find that, 
assuming a Salpeter IMF, the effect of SN explosions is not very 
important, so that omission in the earlier is justified. The main con- 
clusions in lRicotti et alJ i2002d) remain unchanged if we include 
SN feedback and the radiative feedback is strong enough {i.e., as- 
suming a Salpeter IMF and (/esc ) ^ 10%). However, in the current 
work assuming a top-heavy IMF, the SN energy injection has more 
dramatic effects and also affects galaxies with masses A/* ^ 10^ 
Mq. 

ficult to compare predictions of different works (numerical and analytical) 
and estimate which are reasonable values to assume for the free parameter 
(/esc). Nevertheless, (/esc) cannot be treate d purely as a fr eely adjustable 
parameter but has cosmological interest {e. g. ■ lRicottl2003al) . 
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Figure 2. Comoving mean free path Acom of H I ionising photons as a 
function of redshift. The size of H II regions before overlap is approximately 
equal to Acom- Here we show Acom for several simulations in table [Tito 
point out the effects of changing the box size and mass resolution. 



3.2 Results 

In table □ we list the parameters adopted in the simulations. Here 
e« is the star form ation efficiency of the adopted star formation law 
iCen & OstrikeillT992h : dp^/dt — e,p^/t,, where p, and pg are 
the stellar and gas density, respectively. The quantity t, is the max- 
imum of the local dynamical and the local cooling time. The escape 
fraction from a cell, (/esc), is resolution-dependent. The parameter 
FiMF is proportional to the mean metallicity yields and energy in- 
put by SN explosions of the stellar populations. This parameter is 
unity for a Salpeter IMF and population II stars. But Jimf > 1 if 
the IMF is top-heavy and pair-instability SNe are dominant. Pop- 
ulation III stars could have Fimf < 1 if super-massive stars, that 
collapse directly into BHs without exploding as SNe, dominate the 
mass function. 



3.2.1 Convergence study: (/c. 
minihaloes 



and photoevaporation of 



We have run a set of simulations with box sizes Lbox = 1,2 and 
4 h~^ Mpc and number of particles A^^ox = 64^ and A^box = 
128'^. The other parameters of the simulations, listed in table Q 
are the same: we use the maximum efficiency of UV production 
allowed by thermonuclear reactions, euv" — 2 x 10~^, and (/esc) 
= 0.5. We include the effect of SN explosions with an energy input 
and metal production 17 times the value for a Salpeter IMF. This 
is a reasonable assumption for a top-heavy IMF, if pair-instability 
SNe are important. But, given the uncertainties in the energy of 
SN explosions and on the IMF of population III stars, a scenario 
is also plausible in which population III stars collapse directly into 
BHs or produce subluminous SN explosions ejecting few metals. 
We explore this second scenario running a simulation with box size 



Figure 3. Clumping factor of the IGM, C (thick lines), and inside Hll 
regions, Chm (thin lines), as a function of redshift. The solid, dashed and 
dotted lines show the clumping factor, C, for the strong SN feedback sim- 
ulations with A'^box = 128 and Lbox = 1,2 and 4 h~^ Mpc. respectively, 
listed in table The long-dashed lines refer to the Lbox = 2 h~^ Mpc 
simulation with weak SN feedback. 
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Figure 4. Optical depth to Thomson scattering, Tc , as a function of red- 
shift for the simulations with A^box = 128 listed in tableQ The horizon- 
tal lines show the best value of Te and la confidence limits from WMAP 
iKogut et all2003l) . The vertical band indicates the redshift probed by the 
furthest quasars found by the SLOAN survey. 



© 2002 RAS, MNRAS OOO.rilfnl 



8 M. Ricotti and J.P. Ostriker 



Lhox ~ 2 (128L2VMb in table0 with reduced energy input and 
metal ejection (Fimf ~ !)■ 

In the top panels of Fig. Q we show the global star formation 
rate (SFR) as a function of redshift. In each panel, from left to right, 
we show simulations with box size ibox = 1,2 and 4 h^^ Mpc. 
The thick lines show the A'box = 128 simulations and the thin 
lines the A^box = 64 simulations. The simulation with weak feed- 
back is shown in the central panel with a dashed thick line. In the 
high-resolution Lbox = 1 h^^ Mpc simulation we resolve star for- 
mation in galaxies w ith masses Mdm ~ 3.94 x 10^ Mq (see 
iRicotti et alj2002'3) . But, because of the strong feedback from SN 
explosions, star formation in the small-mass galaxies is suppressed. 
This can be seen comparing the top left panel to the simulations 
with larger box sizes and lower mass resolution. In the larger box 
simulations the first stars form later but the SFR quickly converges 
to the same value. The SFRs in the Lbox = 2 and Lbox = 4 h^^ 
Mpc box are almost identical. At z ~ 7 the global SFR is 3 x 10"'^ 
M0yr~^ h~^ Mpc This value is still very small when compared 
to SFR~ 0.1 Mi7)Vr~^ h~^ Mpc estimated from observations 
at z ~ 4 - 5 iLanzetta et all2002h . 

In the bottom panels of Fig. Q we show the integrated num- 
ber of ionising photons per free electron, Nph/ (xc), as a function 
of redshift. The redshifts after reionisation (i.e., (xc) > 0.9) are 
marked with points. The simulations shown are the same as in the 
top panels. At the redshift of reionisation the quantity Ai'ph mea- 
sures the number of recombinations per baryon. If recombinations 
are negligible, A'ph is unity. From the plots it appears that the re- 
combinations in the A'box ~ 64 simulations (thin lines) are always 
underestimated with respect to the simulations with higher mass 
resolution A^box ~ 128 (thick lines). We see that in the smaller 
scale boxes, which have higher mass resolution, and hence allow 
for lower mass halo formation, the integrated number of ionising 
photons produced per free electron required to reionise the IGM 
is larger. This number initially increases as the Strijmgren spheres 
are finding their way through the denser regions surrounding the 
ionising sources. When the volume filling factor of the Hll re- 
gions is larger, the ionisation fronts propagate mostly in the under- 
dense regions and the integrated mean number of recombinations 
per baryon, Nph/{xc), decreases to less than 10. In the left panel 
the H II regions never break out of the dense filaments, because of 
the bursting mode of star formation in small-mass galaxies and the 
absence of large ones in the small box. Because of the self-limiting 
nature of the stellar feedback, as already found by Ricotti et al. 
j2002bh . small-mass galaxies are not able to reionise the IGM. By 
comparing the middle and right panels we see that the integrated 
mean number of recombinations per baryon at the redshift of reion- 
isation (marked by the first dot) is lower for the larger box size, that 
has lower mass resolution. This number is the same for simulations 
that have the same mass resolution (e.g., the thin line in the middle 
panel and the thick line in the right panel). 

We have shown that the integrated number of recombinations 
at the redshift of reionisation depends on the mass resolution of the 
simulation. This means that low resolution simulations might un- 
derestimate the clumping factor inside the H II regions (cf, Fig.|3j 
or overestimate the mean value of (/esc). Resolving small-mass 
haloes, that are photoevaporated during cosmic reionisation con- 
suming a fraction of the ionising photon s, may a lso be impor- 
tant to achieve a c onvergent value of Tc iShapiro & Raga 20o9 ; 
iHaiman et ai]|200ll: lBarkana & Loeb 200^). We also note that the 
bursting mode of star formation in small-mass galaxies increases 
the number of recombinations and thus increases the number of 
ionising photons needed to reionise the IGM. 
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Figure 5. Hydrogen ionisation fraction as a function of redshift for the 
simulations with A'^box = 128 listed in table The circles, squares and 
crosses show the redshift at which the mass weighted mean metallicity in 



each simulation is Z = 10 



, 10- 



, and 10 Zq , respectively. 



In Fig.|2|we show the comoving mean free path Acom of H I 
ionising photons as a function of redshift for the simulations in ta- 
ble Q with strong feedback from SN explosions. The simulations 
with Lbox = 2 h~^ Mpc and Lbox = 4 h"^ Mpc have approx- 
imately the same SFR but in the larger box reionisation (defined 
as the redshift at which Acom increases steeply) happens much ear- 
lier: Zici ~ 7.5 in the L = 2 h^^ Mpc run and z^ci ~ 9.5 in the 
Lbox = 4 h"^ Mpc run. The two Lbox = 2 h~^ Mpc simulations 
with different mass resolutions (A^box = 64 and 128) reionise the 
IGM at the same redshift even if the SFR is half in the smaller 
resolution simulation. This effect can also be seen comparing the 
Lbox = 1 h^^ Mpc simulations with the others. These results show 
that, even if we do not need to resolve small-mass haloes to have 
a convergent global SFR, we do need a high mass resolution or 
otherwise the optical depth to electron Thomson scattering will be 
overestimated due to neglect of recombination in high density re- 
gions. 



3.2.2 Strong and weak SN feedback 

In Fig.|4|we show the optical depth to Thomson scattering. To, as a 
function of redshift for the simulations with A'^box ~ 128 in tableQ 
In Fig.|5|we show the hydrogen ionisation fraction as a function of 
redshift for the same simulations. If we define z^d as the redshift at 
which the hydrogen ionisation fraction is 90%, we find z^ci ~ 7.5 
for the L — 2 h^^ Mpc simulation with strong SN feedback and 
Zici ~ 10 for the simulation with weak SN feedback. For the same 
simulations the optical depths are Tc = 0.09 and Tc = 0.13, respec- 
tively. These optical depths are considerably larger than expected 
for a sudden reionisation at Zrci. This is because the duration of 
reionisation, from the redshift where the mass-weighted mean hy- 
drogen fraction is 10 %, to complete reionisation, is quite long. We 
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Figure 6. (left) Mass weighted metallicity-overdensity distribution for tlie simulation 128L1VM at redshift 2 = 9. The presence of metal-enriched wind with 
Z ~ 10~^Zq, produced by SN explosions in small-mass galaxies (M > 4 X 10® M©), is evident in both the metallicity distribution and density distribution, 
(right) The same as in the left panel but for the simulation 128L4VM. The galaxies resolved in this simulation are less affected by SN explosions because of 
their larger masses (M > 2.5 X 10* Mq). 



find approximately Azrci/^^rci 50 %. The maximum Tc 0.13 
found for the Lbox = 2 h^^ Mpc simulation is obtained for the 
case when SN feedback is not important. But this condition is not 
sufficient. The maximum optical depth, Tc ~ 0.13, can be achieved 
only if the metallicity of the gas in star forming regions remains 
lower than a critical value. In the next section we show that this 
requirement puts strong constraints on the yield that population III 
stars may have in order to contribute to the IGM reionisation. 

In summary, we do not achieve convergence for the simula- 
tions with strong feedback. Rare massive galaxies in this case con- 
tribute to the global SFR, because star formation in the smaller ones 
is suppressed by SN explosions. For this reason we need a box size 
of at least Lbox = 4 h~^ Mpc to achieve a converged estimate for 
the SFR. But when we increase the box size the mass resolution 
becomes too small and we underestimate the clumping of the IGM 
(i.e., the number of recombinations). For this reason, even if the 
SFR has converged, Tc is likely to be overestimated because of the 
too low recombination rate in the Lbox = 4 h^^ Mpc box. The 
ibox = 2 h~^ Mpc run instead gives a better estimate of Tc be- 
cause the underestimated clumping of the IGM is compensated by 
the slightly underestimated SFR. This can be seen by comparing 
the two I/box = 2 h^^ Mpc runs with 64^ and 128"^ particles in the 
middle panels of Fig.Q in the low resolution run the number of re- 
combinations (bottom panel) and the SFR rate are underestimated 
(top panel). Fortunately these two effects almost cancel each other 
and, even if the low resolution simulation (thin lines) is not conver- 
gent, we did get an estimate of the redshift of reionisation (marked 
by the last dot in the bottom panel) and Tc that is about the same as 
in the higher resolution simulation (thick lines). 

For the weak feedback case we should be close to convergence 
since we can safely use the smaller box. In this case the SFR is 
large enough that the contribution from rare massive galaxies is 



negligible. Unresolved small-mass galaxies do not contribute sub- 
stantially to the SFR. Ideally, we need larger mass resolution to 
be sure that we are not underestimating recombinations. Note that 
(/esc) is defined as the fraction of ionising photons that escapes 
from the resolution element. Therefore, by definition, it is resolu- 
tion dependent and as we increase the resolution of the simulations 
its value should approach unity. We would use (/esc) ~ 1 if we 
could perfectly resolve galaxies down to the star forming regions 
in the ISM. Note that by definition, the effective clumping factor of 
the IGM is complementary to (/esc), in the sense that it accounts 
for all the recombinations produced in regions larger or equal to the 
resolution element. 



3.2.3 Metal enrichment 

The circles, squares and crosses in Fig.|5|show the redshift at which 
the mass weighted mean metallicity in each simulation is Z — 
10~^, 10~^, and lO"** Zq, respectively. The volume weighted 
mean metallicities (not shown here) are ten times smaller than the 
mass weighted mean for both the simulations with strong and weak 
SN feedback. This means that most metals are located in overdense 
regions. This can be better seen in Fig.|6| where we show the mass 
weighted metallicity distribution as a function of the overdensity 
for the simulation 128L1VM with Lbox = 1 h~^ Mpc (left) and 
128L4VM Lbox = 4 h"^ Mpc (right). In both simulations the 
mean metallicity increases with the overdensity. There is some scat- 
ter around the mean but at overdensities larger than 200 most of the 
gas has metallicity Z ~ 0.1 — 0.01 Zq (cf., Fig.0. These values 
are about ten times larger than the mass-weighted mean metallici- 
ties at the same redshift (or a hundred times the volume-weighted 
mean metallicities). Similar results were found by Cen & Ostriiffi^ 
("199^. If Zci ~ 10~* - 10~^ is the critical metallicity at which 
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Table 2. Results of simulations with radiative transfer. 
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Figure 7. Metallicity distribution of the gas mass fraction with overdensity 
larger than 200 at redshift 2 = 9. The solid line refers to the simulation 
with box size Lbox = 1 h" ^ Mpc and the dashed hne to the simulation 
with box size Lbox = 4 h~^ Mpc. 



of their largermasses (Mdm > 2.5 x 10 M©), are less affected by 
SN explosions and feedback. 

In table |2| we give the values of relevant quantities measured 
in our high resolution simulations with parameters listed in tableQ 
Here a;, is the fraction of baryons in stars, Zrd is the Gunn-Peterson 
reionisation redshift and A'^ph is the number of ionising photons per 
baryon emitted at the redshift of reionisation. We show the mean 
mass weighted and volume weighted metallicities. In the last col- 
umn we list the mean mass weighted metallicity in regions with 
overdensities > 200. 

As noted in § Q the simple estimate presented in that section 
for the number of ionising photons TVph as a function of the optical 
depth Tc, 



ph 



10 ( — 

Vo.i 



(16) 



with a = 1 (cf., equation (6|), is not very accurate when compared 
to the simulations. A better agreement with the simulations can be 
obtained using instead a = 4. The reason for this result is that 
the mean effective clumping factor inside the H II regions it is not 
constant with redshift as was initially assumed. Using a — 4 and 
g = 2 in equation <13> . the estimated IGM metallicity as a func- 
tion of To is in good agreement with the simulation results listed in 
table|2| Since in equation <13> we have assumed that all the metals 
are ejected from the galactic haloes we should compare To (given 
in the third column) to the mass weighted metallicities (given in the 
seventh column). 



population II stars dominate over population III stars, it is quite 
clear from Fig. |5| that reionisation by population III stars cannot 
be completed in any of our simulations. Since stars are formed in 
overdense regions where the metallicities are larger, population III 
stars can not be important for reionisation unless their yield is much 
smaller than assumed here or the metals do not mix efficiently with 
metal-free gas. Note that we require a factor of 10^ reduction in 
either yield or mixing efficiency to keep the metallicity below the 
level required for production of population III stars (cf.. Fig. 0. 
While Fig.Qimplies that population III is probably not relevant to 
reionisation it is not intrinsically implausible or inconsistent with 
observations. The quasar broad line emission regions sample the 
high density gas in galaxies and the Sloan high redshift quasars 
(z ~ 6) already show a metallicity approaching solar values. This 
is consistent with the most relevant (Lbox = 4 h^^ Mpc) curve in 
Fig.Qand would remain so even if we allow for 50 % fall-back of 
metals onto black holes. 

In the simulation with Lbox = 1 h"^ Mpc, the presence of 
metal-enriched galactic winds with Z ~ lO"* Z©, produced by SN 
explosions in small-mass galaxies (M > 4x10" M©), is evident in 
both the perturbed metallicity distribution and density distribution 
(cf., Fig.|6j. In the Lbox = 4 h"^ Mpc simulation galactic outflows 
are not important. The galaxies resolved in this simulation, because 



4 SEMIANALYTIC SIMULATIONS 

In this section we use a simple semianalytic model to study the 
reionisation history of the IGM. This model can be useful to study 
the dependence of reionisation history on the assumed cosmologi- 
cal model and parameters. The semianalytic model do not include 
the effects of radiative feedback, SN feedbacks and metal enrich- 
ment. It also uses the results of the numerical simulations on the 
clumping factor of the IGM. Given these assumptions and limita- 
tions (due to neglecting feedbacks) the results of the semianalytic 
models are in good agreement with the numerical simulations. For 
example, the value of Tc — 0.13 found in our numerical simulation 
with weak SN feedback (128L1 VMb) is shown as a circle in Fig.|8| 
The corresponding value found in the semianalytic model, shown 
with a square, is about 15% larger. 

4.1 The code 

We implemented a semianalytic model to study reionisation, chem- 
ical evolution and re-heating of the IGM. The method to calcu- 
late the filling fact or of H ll regions and the SFR is the same as 
in Chiu & O striked ilOOOl) . But we also include X-ray ionisation of 
the IGM before complete overlap of the H II regions. Since in this 
paper we neglect the X-ray background, a full description of the 
code will be given in Appendix A of paper II. Here we only briefly 
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Figure 8. Optical depth to Thompson scattering, Tc, as a function of the 
spectral index of initial density perturbations, n^, computed from the semi- 
analytic code. Here erg is calculated according to the amplitude of per- 
turbations at fc = 0.05 h Mpc~^ measured by WMAP. Each line from 
top to bottom shows Te iris) for effective efficiencies of UV emission 
= 10~^, 10 ~^ an d 10~*, respectively. We show the Icr confi dence 
contour of Te-Us iSpe rgel et al. 2003) and of Tc ' Koaut et al. 2 0031) . The 
value of Te = 0.13 found in our numerical simulation with weak SN feed- 
back (128LlVMb) is shown as a circle. The corresponding value found in 
the semianalytic model, shown with a square, is about 15% larger. 



summarise the main features of the semianalytic model. The mass 
function of DM haloes and their formation/merger rates are cal- 
culated using the extended Press-Schechter formalism. The SFR 
is assumed to be proportional to the formation rate of haloes. The 
effects of cooling and dynamical biases that depend on the mass 
of the haloes are taken into account. We consider the IGM as a 
two-phase medium: one phase is the ionised gas inside the H II re- 
gions and the other is the neutral or partially ionised gas outside. 
We solve the radiative transfer for the volume averaged specific 
intensity and we derive the specific intensity inside and outside the 
H II regions from their volume filling factor and separating the con- 
tribution from local UV sources and the background radiation from 
distant sources. We solve the time-dependent chemical network for 
eight ions (H, H+, H", H2, H^, He, He+, He++), and the ther- 
mal evolution as a function of the gas density outside the H II re- 
gions. Given the clumping factor of fully ionised gas around the 
UV sources (taken from the SLH simulations) we evolve the filling 
factor, temperature and chemistry inside the Hll regions. The fol- 
lowing heating and cooling processes are included: collisional- and 
photo-ionisation heating, ionisation by secondary electrons, H, He, 
H2 cooling, and Compton and adiabatic cosmological cooling. The 
rates are the same as in cRicotti et al. 2001) . 



4.2 Dependence of Tc on small-scales power 

In this section we use the semian alytic code to study the depen- 
dence of Te on the cosmology fe.g.. lYoshida et all2003bl£l) . In par- 
ticular we focus on models with a power spectrum of perturbations 
that differs in power at small-scales. Here we study the dependence 
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Figure 9. Optical depth to Thomson scattering, Te, as a function of the 
effective efficiency of UV emission (/esc) ^uv f™ CDM (solid) and a 
WDM (dashed) with particle mass 1.25 MeV cosmologies. We use a semi- 
analytic model with free parameters chosen to reproduce the results of the 
numerical simulation and the observed SFR at z ^ 6. Here feedback effects 
are neglected. We estimate that they would reduce Te by 20-25%, as shown 
in the shaded region for the CDM model, and leave the WDM model unal- 
tered. The vertical solid line shows the maximum effective efficiency avail- 
able from thermonuclear sources. Plausible (population II) models would 
shift it to the dotted line reducing To to the level of 7-8% which is the value 
directly inferred by observations of quasar reionisation <Fan et alj Eool 
IChiu et all2003h and is below the 2cr error bar of WMAP. 



of Te on the slope of the initial power spectrum of perturbations 
for a CDM and warm dark matter (WDM) cosmology. The results 
presented in this section are intended to show the dependence of 
Tc on the cosmological parameters, but they are not sophisticated 
enough to give a precise estimate of Tc. The reader should keep in 
mind that in these models the most important feedback processes - 
those due to SNe - are neglected. In Fig.|8|we show Te as a function 
of the slope, n^, of the power spectrum of initial density pertur- 
bations. We normalise the power spectrum according to the ampli- 
tude A at wave nu mber k = 0.05 h Mpc~^ measured by WMAP 
jVerde et alj|2003h . Consequently the variance ag at z = is de- 
termined by the value of Us. We show the results for values of the 
product (/csc) euv = = 10"^ 10"^ and 10^". The shaded 
contours show the Icr confidence limits in the Tc-ris plane derived 
by the WMAP ana l ysis of temperatur e fluctuations of the CMB 
JSoergel et alj2003tKogut et alfcOOsI) . A small tilt, Us < I, sup- 
presses the power on small scales delaying star formatio n in the first 
galaxies and the redshift of reionisation, as also found bv lChiu et alJ 
(.2001) . We find a linear relation between Te and the spectral index 



Te ~Te(n, = l)4.54(n, -0.78), 

where Te{ns — 1) is the value of Tc for Us — 1. In Fig.|9|we 
show Te (ris = 1) as a function of ey'v for the ACDM and AWDM 
cosmology. The optical depth Tc increases exponentially with euy: 



Teiris 



a + b log 



eff 
10-3 
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where a = 0.145(0.10) and h = 0.055(0.036) for the ACDM 
(A WDM ) cosmology. Our results roughly agree with a similar 
study by IChiu et all i2003h . but a precise comparison with their 
work is not possible since we adopt a constant effective efficiency 
of UV emission, ejjv, while they adopt a time dependent effective 
efficiency. 

The value of Tc = 0. 13 found in our numerical simulation with 
weak SN feedback (128LlVMb) is shown as a circle in Fig.[8| The 
corresponding value found in the semianalytic model, shown with a 
square, is about 15% larger. This can be due to the finite resolution 
of the numerical simulation or the radiative and weak SN explo- 
sion feedbacks included in the numerical simulation but, consider- 
ing these difference in the underlying assumptions, the agreement 
between the two methods is encouraging. In the numerical simula- 
tions including strong feedback, the SN explosion heating reduces 
the value of Tc by 20 — 25% as compared to the weak feedback 
case. 

We have also examined a warm dark matter model with our 
semianalytic set of calculations. If the strong SN feedback seen in 
the numerical simulations produces a similar reduction of To in the 
semianalytic CDM model, we find that the warm dark matter model 
with particle mass 1.25 keV does, as expected, produce a lower re 
but the effect is surprisingly small and estimated to be on 10 % 
for a given value of (cf., Fig.|9}. The effect of SN explosions 
will leave the WDM model essentially unchanged because small- 
mass galaxies are not present in this model. The photoevaporation 
of minihaloes is expected to reduce Tc in CDM models, but this 
effect is not present in WDM models. 



5 CONCLUSION 

We have used cosmological hydro-dynamical simulations with ra- 
diative and SN explosion feedback to study the reionisation of the 
IGM by population m stars. Simple calculations have been used 
to argue that a large Tc is difficult to achieve and is incompatible 
with low metal enrichment of the IGM/ISM if population III stars 
end their lives as pair-instability SNe or smaller-mass hypemovae. 
A semianalytic code is used to study the dependence of Te on cos- 
mological parameters. 

The main results from the cosmological simulations are sum- 
marised in the following points. 

(i) It is very difficult or impossible to use stellar sources to 
reionise the universe at Te > 0.13 given a maximum efficiency 
of ionising radiation production typical (extreme population III) 
of thermonuclear fusion reactions eyv^ = 2 X 10-^ and (/esc) 
= 0.5. When feedback from SN explosions is included, the lim- 
ited mass resolution of the cosmological simulations is not crucial 
for achieving convergence for the global star formation history in 
the first galaxies. However a larger box size is required because the 
galaxies that contribute mostly to star formation are those that are 
the most massive. Our studies of convergence indicate that increas- 
ing the mass resolution over what we utilised in the I28L2VM run 
(which gave Tc = 0.13) would not increase Tc for models with SN 
feedback. 

(ii) If metal production is normal (top-heavy IMF, or pair in- 
stability SNe or hypemovae), the mechanical energy input by SN 
explosions is about ten times larger than for Salpeter IMF. This pro- 
duces strong outflows in galaxies with masses Mdm ;s io» M© , re- 
ducing their star formation and delaying reionisation to Zrci < 10. 

(iii) If the ratio of metal to ionising photons from metal poor 
population III stars is normal, then metal enrichment of the ISM 



prevents metal-poor stars from being produced for long enough to 
reionise the IGM. One requires a factor of 10'^ reduction of yield 
or efficiency of mixing to alter this conclusion. 

(iv) The issues raised in the previous two points can be allevi- 
ated if most population III stars collapse into a BH without pro- 
ducing a large amount of metals. This solves the metal pollution 
problem making the epoch of population III star domination longer. 
And in this scenario the mechanical feedback from SN/hypernovae 
explosions would also be reduced. 

(v) If metal poor stars are initially important, then the secondary 
production of ionising radiation due to BH accretion may be more 
important than the primary production. 

(vi) As a byproduct of the semianalytic treatment we find that 
a warm dark matter model with particle mass 1.25 keV does, as 
expected, produce a lower Tc, but the effect is surprisingly small, 
and estimated to be on 10 % for a given value of tyv. 

(vii) If we (self consistently) neglect population III and have 
more normal population II star formation properties we derive a 
much lower Tc of 0.07 — 0.08 consistent with other authors but 
inconsistent with WMAP results. 

In conclusion, if the first stars were massive and most of the 
radiation escaped from the host galaxies, the increase of Tc could 
be within the la confidence limits of WMAP data. But these two 
conditions are not sufficient. The energy liberated by SN explo- 
sions and the metal yields would have to have been much smaller 
than expected. This requirement seems to agree nicely with one 
possible interpretation of the abundance pattern in the recently dis- 
covered most iron-deficient star, HEO 1 07-5240, and other metal- 
poor but carbon-rich stars. Subluminous supernova explosions with 
E51 ~ 0.3 are characterised by a large fall-back on the central 
black hole and their ejecta are carbon rich and iron deficient. The 
gas, enriched in carbon and oxygen, can cool fast enough to pro- 
duce a second generation of star s with abu ndance patterns similar 
to HE0107-5240 jumeda & No motdl20o3) . In addition the abun- 
dance pattern of HE0I07-5240 is not compatible with the yields 
from pair-instability SNe. An interesting suggestion is that the en- 
ergy of a SN might be low if the BH remnant is not spinning, 
and large for spinning BHs. Note that these conclusions are true 
even if the first stars were not supermassive. Indeed, the abundance 
pattern of HE0107-5240 can be explained assuming a mass of 25 
M0 for the SN progenitor which is roughly the mass indicated by 
I Abel et ali 12002) calculation of "the first star". An interesting con- 
sequence of this scenario is a copious production of rather massive 
BHs from the first stars. For this reason, the secondary radiation 
from accretion on seed BHs might be a more important source of 
ionising radiation than the primary radiation from the massive stars 
(see also. lMadau et al . 2003). BH accretion quickly builds up an X- 
ray background that can keep the IGM partially ionised before the 
complete reionisation by stellar sources at z ~ 7. We investigate 
this scenario in separate papers (paper Ila and paper lib). 
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